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TECHNICAL  NOTE  29l6 

EFFECT  OF  THERMAL  PROPERTIES  ON 
LAMINAR-BOUNDARY-LAYER  CHARACTERISTICS 
By  E.  B.  Klunker.  and  F.  Edward  McLean 

SUMMARY 


An  iteration  method  is  presented  for  solving  the  laminar -boundary- 
layer  equations  for  compressible  flow  in  the  absence  of  a  pressure 
gradient  wherein  the  temperature  variation  of  all  the  fluid  thermal 
properties  is  considered.  Friction  and  heat-transfer  characteristics 
have  been  calculated  for  a  stream  temperature  of  -67°  F  for  Mach  numbers 
from  1  to  10  with  the  use  of  values  of  the  heat  capacity,  conductivity, 
and  viscosity  determined  from  experiment.  Consideration  of  the  temper¬ 
ature  variation  of  all  the  fluid  thermal  properties  causes  the  recovery 
factor  to  decrease  substantially  with  increasing  Mach  number.  Moreover, 
the  heat-transfer  rate  is  found  to  be  proportional  to  the  difference 
between  an  effective  enthalpy,  which  is  a  function  of  both  the  surface 
temperature  and  stream  Mach  number,  and  the  surface  enthalpy.  In  con¬ 
trast,  the  heat-transfer  rate  is  approximately  proportional  to  the  dif¬ 
ference  between  the  recovery  enthalpy  and  the  surface  enthalpy  for 
solutions  which  employ  a  constant  Prandtl  number.  The  calculated  skin 
friction  and  heat-transfer  rates  based  upon  the  use  of  the  Sutherland 
equation  for  viscosity  and  a  Prandtl  number  of  0. 75 >  however,  are  in 
excellent  agreement  with  the  results  of  the  present  analysis. 


INTRODUCTION 


Most  compressible -flow  boundary- layer  analyses  employ  approximate 
relations  for  the  temperature  variation  of  the  thermal  properties .  A 
simple  analytic  relation  is  usually  specified  for  the  variation  of  vis¬ 
cosity  with  temperature.  Moreover,  the  heat  capacity  and  the  Prandtl 
number  are  generally  considered  to  be  independent  of  temperature  and 
as  a  consequence  the  conductivity  and  viscosity  are  required  to  have 
the  same  temperature  variation.  Experiment  shows  that  these  conditions 
are  not  satisfied  closely  over  the  wide  temperature  range  which  occurs 
in  a  boundary  layer  in  a  high-speed  flow.  The  effect  of  these  approxi¬ 
mations  on  the  calculated  boundary-layer  characteristics  is  most  easily 
studied  by  considering  the  laminar -boundary-layer  flow  on  a  flat  plate. 
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Calculations  of  boundary-layer  characteristics  which  employ  approxi¬ 
mate  relations  for  the  thermal  properties  are  generally  satisfactory  as 
long  as  the  temperature  variation  in  the  boundary  layer  is  not  excessive; 
the  range  for  which  the  solutions  are  satisfactory  depends  upon  the  value 
of  the  Prandtl  number  used  in  the  analysis  and  upon  the  viscosity- 
temperature  relation  employed.  Various  viscosity-temperature  relations 
and  values  of  the  Prandtl  number  have  been  used  in  the  analysis  of  the 
laminar  boundary  layer  on  a  flat  plate  (refs.  1  to  5,  for  example) .  The 
choice  of  a  Prandtl  number  of  unity  (ref.  1,  for  example)  or  of  a  linear 
viscosity -temperature  relation  (ref.  2,  for  example)  leads  to  an  essen¬ 
tial  simplification,  since  in  each  case  only  one  differential  equation 
need  be  solved.  Better  solutions  are  obtained  by  taking  the  Prandtl 
number  as  approximately  0 .73  (the  value  for  air  at  normal  temperatures) 
and  by  utilizing  a  more  realistic  viscosity- temperature  relationship.. 

The  viscosity-temperature  relationship  is  frequently  taken  as  a  simple 
power  law  (refs.  1  and  3 ,  for  example)  or  as  given  by  the  Sutherland 
equation  (refs.  4  and  5,  for  example) .  The  Sutherland  equation  for  vis¬ 
cosity  shows  especially  good  agreement  with  experiment  over  a  wide  temper¬ 
ature  range. 

At  high  supersonic  Mach  numbers,  where  the  temperature  variation 
within  the  boundary  layer  is  large,  the  assumption  of  a  constant  Prandtl 
number  and  the  use  of  an  approximate  viscosity-temperature  relation  is 
not  jusitif ed  a  priori .  A  few  calculations  of  the  laminar -boundary- 
layer  characteristics  on  a  flat  plate  have  been  made  which  have  not 
imposed  these  restrictions.  Experimental  values  of  the  thermal  properties 
were  utilized  in  some  calculations  of  the  laminar -boundary- layer  charac¬ 
teristics  by  an  iteration  method  at  Mach  numbers  from  1  to  10  in  refer¬ 
ence  6  and  the  thermal  properties  employed  in  the  differential-analyzer 
solutions  of  references  7  and  8  were  taken  in  part  from  experiment  -  For 
the  latter  two  analyses,  the  stream  conditions  were  determined  from  the 
solution.  Consequently,  a  graphical  or  trial -and-error  procedure  was 
required. in  the  course  of  the  solution  to  obtain  values  corresponding 
approximately  to  the  desired  stream  conditions.  In  the  iteration  methods 
used  in  reference  6  and  in  the  present  analysis,  the  solutions  are  deter¬ 
mined  for  specified  stream  conditions . 

The  present  investigation  is  a  continuation  of  the  work  reported  in 
reference  6.  The  purpose  is  to  determine  accurately  some  characteristics 
of  the  boundary  layer  at  high  supersonic  speeds  and  to  show  the  effect 
of  various  approximations  to  the  thermal  data  of  the  fluid  on  the  accuracy 
of  the  calculated  friction  and  heat-transfer  characteristics.  Tabulated 
values  of  the  thermal  properties,  found  from  experiment,  have  been 
employed  in  the  calculations  in  order  to  obtain  reliable  data.  Friction 
and  heat-transfer  characteristics,  calculated  for  a  stream  temperature 
of  -67°  F  for  Mach  numbers  from  1  to  10,  are  presented  in  the  form  of 
graphs.  Tabulated  values  of  the  corresponding  velocity  and  temperature 
profiles  and  their  derivatives  are  available  upon  request  from  the 
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National  Advisory  Committee  for  Aeronautics.  The  iteration  method 
described  herein  for  the  solution  of  the  equations  is  similar  to  that 
presented  in  reference  6.  The  computations  are  facilitated  however 
through  a  different  formulation  of  the  equations;  the  convergence  of 
the  iteration  procedure  is  improved  considerably  and  a  troublesome 
singularity  is  avoided. 

SYMBOLS 


x>y 

u,v 

P 

V- 

i 

k 

T 


cp 


v 


ui 

Pr 


coordinates  parallel  and  normal  to  stream  direction 
components  of  velocity  along  x-  and  y-axes 
density 

coefficient  of  viscosity 

enthalpy  per  unit  mass 

thermal  conductivity 

absolute  temper ature 

heat  capacity  at  constant  pressure 

kinematic  viscosity 
stream  velocity 

Prandtl  number 


$ 

<X,P 


enthalpy  function. 


integrating  factors 


k 
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T 

R 


shear  stress  at  surface 
Reynolds  number 

average  skin-friction  coefficient 


1 

7 

M 

r 

c 

Subscripts : 


local  heat-transfer  rate  per  unit  area  at  surface 

ratio  of  heat  capacities  at  constant  pressure  and  con¬ 
stant  volume 

Mach  number 

temperature  recovery  factor 
constant 


free  stream  or  at  edge  of  boundary  layer 
stagnation  condition  of  stream 
effective 
r  recovery 

s  surface 

Superscript: 

*  dimensionless  quantity  based  on  stream  conditions 


1 

0 

e 


ANALYSIS 


Basic  equations,-  With  the  x-axis  taken  in  the  free-stream  direc¬ 
tion,,  the  two-dimensional  compressible-flow  boundary-layer  equations 
for  steady  motion  in  the  absence  of  a  pressure  gradient  may  be  written 
as  follows : 


Momentum  equation: 
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Continuity  equation: 


Energy  equation: 


^-(pu)  +  ^-(pv)  =  0 
ox  oy 


pu  — +  pv 

dx 


i  ,  9-4  + 

y  9y\  9y/ 


where  u  and.  v  are  the  velocity  components  in  the  direction  of  the 
x-  and  y-axes,  respectively,  p  is  the  density,  T  is  the  absolute 
temperature,  i  is  the  enthalpy  per  unit  mass,  H  is  the  coefficient 
of  viscosity,  and  k  is  the  thermal  conductivity.  The  heat -conduct ion 
term  of  equation  (3)  may  be  expressed  in  terms  of  the  enthalpy.  With 

the  use  of  the  relation  c„  =  where  c_  is  the  heat  capacity  at 

y  ST  * 

constant  pressure,  the  energy  equation  becomes 


pu  |i 
ox 


JL/JL 

dy  \cp 


y /  V^y  / 


For  a  flow  with  heat  transfer  (specified  constant  surface  temper¬ 
ature)  or  for  a  flow  with  zero  heat  transfer,  the  three  partial- 
differential  equations  which  describe  the  motion  in  the  boundary  layer 
may  be  reduced  to  a  set  of  two  ordinary  differential  equations . '  For 


this  purpose  the  laminar -boundary-layer  parameter  T)  =  y 


where 


Uq  is  the  stream  velocity  and  Vq  is  the  kinematic  viscosity  based 

on  the  stream  temperature,  is  taken  as  the  independent  variable.  Trans 
forming  equation  (2)  and  solving  for  pv  results  in 


1 

pv  =  r 


T)  rr-(  pu)  dTj 

dii 
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and  after  an  integration  by  parts 


(5) 


Transforming  equations  (l)  and  (1+)  and  making 


use  of  equation  (5)  gives 


vldT|\ 


Pu  dnj  ji  + 


=  0 


With  the  introduction  of  the  Prandtl  number  Pr,  the  nondimens ional 
function 


■6  = 


1  -  ii  1  -  *1 


*0  "  al  Ux2/2 


and  the  dimensionless  quantities 


*  u 
u  =  — 
U*i 


pi 


*  M- 


Pr  =  — 
Pr- 


equations  (6)  and  (7)  become 
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« 


and 


ffe  ^  f)  - 

dT]\Pr  dT]/  M-  \Pr  dtj  /  \cLt)  / 


(9) 


where 


f (ti)  =  -  f  pV  dTj 

2  J  0 

The  subscript  1  refers  to  values  at  the  edge  of  the  boundary  layer 
or  values  in  the  free  stream  and  ±q  is  the  stagnation  enthalpy  of 

the  stream. 

The  transformation  of  the  boundary- layer  equations  employed  herein 
follows  essentially  that  of  Schuh  (ref.  9)*  This  formulation  of  the 
boundary-layer  equations  for  the  flat  plate  has  certain  advantages  over 
those  in  which  the  velocity  (ref s .  4  and  5,  for  example)  or  a  function 
related  to  the  stream  function  (ref.  6,  for  example)  is  used  as  the 
independent  variable.  A  singularity  arises  in  the  solution  of  the 
equations  in  the  latter  two  forms  which  complicates  the  numerical  solu¬ 
tion  of  the  equations.  In  contrast,  no  singularity  occurs  in  the  solu¬ 
tion  of  equations  (8)  and  (9) •  Moreover,  the  iteration  method  used  to 
solve  these  equations  converges  satisfactorily. 

Equations  (8)  and  (9)  can  be  solved  numerically  by  a  method  of 
successive  approximations.  The  thermal  quantities  are  functions  only 
of  the  temperature  and  consequently  functions  of  the  independent  vari¬ 
able  T]  •  Therefore,  the  coefficients  of  the  derivatives  in  equations  (8) 
and  (9)  may  be  considered  as  known  functions  of  -  T)  if  some  initial 
solution  is  given.  With  the  coefficients  considered  as  known  functions 
of  the  independent  variable  and  with  the  integrating  factors 


and 
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the  equation  of  motion  (8)  and  the  energy  equation  (9)  may  he  written  as 


Solution  of  boundary- value  problem.*-  The  solution  of  equations  (10) 
and  (ll)  gives  the  distribution  of  velocity  and  enthalpy  throughout  the 
boundary  layer  subject  to  certain  conditions  on  the  boundaries.  Solu¬ 
tions  which  have  zero  heat  transfer  at  the  surface  and  solutions  with 
heat  transfer  (specified  constant  surface  temperature)  are  considered 
herein.  The  boundary  conditions  to  be  satisfied  are: 

For  an  insulated  surface : 

At  the  surface  (tj  =  0),  u*  =  0  and 


At  infinity  (tj  =  00)  ,  u*  =  1  and  $  =  0 
For  a  specified  constant  surface  temperature: 

At  the  surface  (q  =0),  u*  =  0  and  *8  ="0S  =  Constant 
At  infinity  ,  (t)  =  °°)  y  u*  =  ,1  and  $  =  0 
The  subscript  s  denotes  the  values  at  the  surface. 

The  first  integral  of  equation  (10)  is 


* 

*  du  ~  /  \ 

\i  -  =  C  a(T}) 

dt]  x 


4 


Since  the  boundary  conditions  on  the  velocity  are  the  same  for  both 
the  boundary  conditions  of  specified  constant  surface  temperature  and 
zero  heat  transfer  at  the  surface,  equations  (12)  and  (13)  are  valid  for 
both  surface  conditions. 

The  temperature  or  enthalpy  distribution  throughout  the  boundary 
layer  is  obtained  by  integrating  equation  (ll) .  The  first  integral  of 
equation  (ll)  satisfying  the  condition  of  zero  heat  transfer  is 


* 


M*  d-a 
Pr*  dtj 


=  -P(Ti)<pCn) 


(14) 
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where 


<p(n) 


=  2Pr1 


|j*  /du*\ 

P\a'n  / 


2 

in 


and  after  a  second  integration  the  enthalpy  function  $  is  found  as 


(15) 


At  the  surface  $  becomes  the  enthalpy  recovery  factor 


The  first  integral  of  equation  (ll)  for  flow  with  heat  transfer 
(specified  constant  surface  temperature)  is 


where 


-if  M  =  c,p 

Pr*  dl  2 


pep 


and  the  function  C2  is  dependent  on  the  surface  temperature,  stream 

temperature,  and  Mach  number  which  is  to  be  determined  from  the  boundary 
conditions.  A  second  integration  gives 
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where 


u(n)  =  r  ^  P* 

J0  |i* 


p1!  n  * 

v(n)  =  /  ^  P  cltj 

J0  n 


From  the  boundary  condition  at  infinity,  C2  is  determined  as 


•3  _  3 

e„  =  -s - -s. 


where 


<  =  U(“) 


With  this  value  of  C. 


M*  dj)  _  ^e  ^s 
Pr*  dr)  V(°°) 


P  -  P<P 


and  the  enthalpy  function  is 


*0i)  =  *„  -  u(n) 


m if*  -  3 ) 

V(oo)  V  e  s) 


Equations  (13)  and  (15)  give  the  distribution  of  velocity  and 
enthalpy  (and  consequently  the  temperature)  throughout  the  boundary 
layer  for  the  insulated  surface  whereas  equations  (13)  and  (17)  give 
the  distribution  of  velocity  and  temperature  corresponding  to  a  specified 
constant  surface  temperature .  The  distributions  of  shear  and  heat  trans¬ 
fer  are  determined  from  equations  (12)  and  (l 6) . 
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Method  of  numerical  solution.-  The  solution  of  the  equations  for 
the  velocity  or  temperature  distributions  is  one  of  successive  approxi¬ 
mations;  consequently ,  some  approximate , solution  for  the  distribution 
Of  velocity  and  temperature  is  required.  The  Blasius  solution  may  be 
taken  for  the  distribution  of  velocity  and  the  parabolic  relation  between 
temperature  and  velocity  for  a  Prandtl  number  of  unity  (ref.  10)  may  be 
used  to  determine  the  corresponding  temperature  distribution.  Fewer 
iteration  steps  are  required,  however,  if  a  more  refined  solution,  refer¬ 
ences  2  and  5,  for  example,  is  taken  for  the  first  approximation.  With 
the  initial  values  determined  by  the  method  of  reference  2,  four  or 
five  iterations  were  required  to  obtain  differences  of  less  than  l/2  per¬ 
cent  between  the  last  two  iterations  (for  both  the  shear  stress  and  the 
heat-transfer  rate,  at  the  surface)  . 

foliowing  procedure  was  found  to  be  most  satisfactory  in  the 
numerical  solution.  With  the  initial  values  of  velocity  and  temperature 
given  for  a  specified  surface  temperature  and  Mach  number  as  a  function 
of  the  independent  variable  r\  and  with  the  variation  of  viscosity  and 
Prandtl  number  with  temperature  prescribed,  the  functions  f (r)),  0,(11) 

and  J(r])  are  determined  by  numerical  integration  (the  trapezoidal  rule 
was  found  to  be  very  satisfactory).  The  velocity  distribution  is  then 
obtained  from  equation  (13) .  By  making  use  of  the  shear  function  given 
by  equation  (12)  and  the  previously  calculated  values  of  f(ri),  the 
unc  ions  p(n),  <P(t)),  U(tj),  and  V(tj)  are  determined  through  numerical 
tii!wvrt0n  Z  the  enthalpy  function  « (rj)  is  determined  fr!S  equa- 
thP  PnS-i  corresponding  temperature  distribution  may  be  found  from 

the  enthalpy  distribution  with  the  use  of  enthalpy  tables.  With  these 

s  f  “4  temperature,  the  step? outlined  above  “e 

J?^^®d,and  f  ^cond  approximation  to  the  velocity  and  temperature 
distribution  is  determined.  This  process  is  continued  until  the  desired 
accuracy  is  obtained.  The  procedure  for  calculating  the  velocity  and 
temperature  distribution  for  the  insulated  surface  is  similar  to^that 
scribed  for  the  condition  of  constant  surface  temperature.  The  enthalpv 
function  in  this  case,  however,  is  found  from  equation  (15) .  W 

Experimental  values  of  the  thermal  properties  or  values  found  from 

eXpre!f  nS  156  used  with  facility  in  the  computa- 

tionp.  The  equations  may  be  simplified  somewhat  for  the  particular  case 
in  which  the  Prandtl  number  is  taken  as  a  constant;  then,  Pr*  =  1  and 

the  integrating  factor  3  reduces  to  aPr.  For  solutions  involving  a 
constant  Prandtl  number,  however,  the  method  of  reference  5  is  preffrable 
since  the  solution  can  be  obtained  with  less  calculation.  P 


KESULT3  AND  DISCUSSION 


Some  friction  and  heat-transfer  characteristics  at  supersonic  Mach 
numbers  have  been  computed  by  the  method  described  in  the  preceding 
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section.  Faired  experimental  values  of  the  viscosity  and  Prandtl  num¬ 
ber,  taken  primarily  from  reference  11,  have  been  used  in  the  computa¬ 
tions.  (The  calculations  were  limited  to  the  temperature  range  in  which 
experimental  data  were  available.)  The  free-stream  temperature  was  taken 
as  392. 7°  R  (-67°  F),  the  value  at  the  isothermal  region  of  the  MCA 
standard  atmosphere,  and  the  corresponding  values  of  the  thermal  proper¬ 
ties  were  taken  as: 

CP1  =  7.718  Btu/(slug) (deg) 

=3*058  x  10 “7  slugs/(ft) (sec) 
k1=  3.227  x  10"6  Btu/(ft) (sec) (deg) 

Pr-L  =  0.73 

The  boundary-layer  characteristics  presented  herein  include  those 
cases  which  were  calculated  in  reference  6;  the  same  thermal  data  were 
used  in  both  analyses.  Some  small  differences  were  found  between  the 
results  of  reference  6  and  those  calculated  by  the  method  presented 
herein.  These  differences  are  due  primarily  to  the  poor  convergence 
of  the  solutions  in  reference  6.  Additional  calculations  have  been 
made  for  those  cases  by  means  of  the  method  of  this  report  in  order  to 
make  the  results  consistent. 

The  calculated  velocity  and  temperature  profiles  for  several  super¬ 
sonic  stream  Mach  numbers  for  the  insulated  surface  are  presented  in 
figure  l(a) ,  and  profiles  for  Ts*  =1,  2,  and  k  are  presented  in  fig¬ 
ures  l(b),  1(c),  and  l(d) ,  respectively .  (Tabulated  values  of  the 
velocity  and  temperature  distributions  and  their  first  derivatives,  to 
four  decimal  places,  together  with  a  chart  of  the  thermal  properties 
used  herein,  can  be  obtained  upon  request  from  the  National  Advisory 
Committee  for  Aeronautics.)  The  corresponding  friction  and  heat-transfer 
characteristics  are  presented  in  figures  2  to  6  together  with  the  results 
of  several  other  analyses  which  make  use  of  various  approximations  for 
the  viscosity  and  Prandtl  number. 

Skin  friction.-  The  shear  stress  at  the  surface  t  is  given  by 


lit- 
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where  R  is  the  Reynolds  number  Vq .  Substituting  from  equation  (13) 

results  in  the  following  equation  for  shear  stress: 

Plul2  1 

T  = - 

ft  j(») 

since  a(0)  =  1.  The  average  skin-friction  coefficient  cf  is  found 
by  integration  as 

X 

Cf  =  - — r  Tdx  =  i  -1-  (18) 

1  Plux2x  ft  J(“) 

2 


The  skin-friction  coefficients  found  from  equation  (l8)  for  the 
insulated  surface  and  for  the  condition  of  heat  transfer  are  presented 
in  figures  2(a)  and  2(h),  respectively,  together  with  the  results  from 
other  analyses  for  a  stream  temperature  of  -67°  F.  All  the  analyses 
show  substantially  the  same  variation  of  skin  friction  with  Mach  number 
for  the  insulated  surface  except  for  the  solution  based  upon  the  linear 
viscosity  relation  li*  =  T*  for  Pr  =  1.0;  for  these  conditions  the 
skin  friction  is  independent  of  the  Mach  number.  The  present  results 
and  those  of  reference  5,  which  are  based  upon  a  constant  Prandtl  number 
of  0.75  and  the  use  of  the  Sutherland  equation  fbr  viscosity,  are  in 
excellent  agreement  for  both  the  insulated  surface  and  for  the  specified 
constant  surface  temperatures .  The  differences  between  the  results  of 
this  paper,  those  of  Moore  (ref.  8),  and  those  of  Young  and  Janssen 
(ref.  7)^  may  possibly  be  attributed  in  part  to  some  inaccuracies  in 
the  differential-analyzer  calculations,  differences  in  the  thermal  data, 
and  a  'small  difference  in  stream  condition  -  the  calculations  from  refer¬ 
ences  7  an&  8  are  based  upon  a  nominal  stream  temperature  of  400°  R. 

^Tn  reference  8,  experimental  values  of  viscosity  are  used  for  tem¬ 
peratures  up  to  3,240°  R  and  for  Prandtl  number,  up  to  2,000°  R;  for 
temperatures  between  3,240°  R  and  the  dissociation  temperature  the  vis- 

!  v  0  *6 

cosity  was  calculated  with  the  use  of  the  power  law  -rr -  =  (  — — ) 

T=32^0  \3240/ 

and  for  temperatures  over  2,000°  R  the  Prandtl  number  was  determined  from 
calculated  values  of  the  thermal  properties .  Values  of  the  Prandtl  num- 
ber  from  various  sources  and  the  Sutherland  equation  for  viscosity  were 
employed  in  the  boundary-layer  calculations  reported  in  reference  7. 
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The  agreement  of  the  results  based  on  a  power  law  and  a  linear  viscosity- 
temperature  relation  with  the  results  of  this  analysis  is  not  so  satis¬ 
factory  as  that  of  the  results  which  make  use  of  the  Sutherland  equation 
for  viscosity. 

It  should  be  noted,  however,  that  the  variation  of  skin-friction 
coefficient  with  Mach  number  is  dependent  on  the  stream  temperature  and 
thus  the  differences  between  the  results  based  on  a  power  law  and  a 
linear  viscosity-temperature  relation  and  the  results  of  this  analysis 
would  depend  on  the  stream  temperature  at  which  a  comparison  was  made. 

With  a  linear  viscosity- temperature  relation  of  the  form 
H*  =  Constant  T*,  the  skin  friction  is  independent  of  Mach  number; 

taking  the  constant  as  unity  results  in  cf\/R  =  1.328  for  all  sur¬ 
face  temperatures  and  stream  Mach  numbers  irrespective  of  the  value ^of 
the  Prandtl  number.  The  linear  viscosity  relation  H  =  Constant  T 
was  utilized  in  the  calculations  of  reference  2  for  a  Prandtl  number 
of  0.72  and  the  constant  was  adjusted  to  give  the  correct  value  of  the 
viscosity  at  the  surface .  This  procedure  yields  a  variation  of  skin 
friction  with  Mach  number  for  the  insulated  surface  since  the  recovery 
temperature  varies  with  Mach  number;  for  the  condition  of  constant  sur¬ 
face  temperature,  the  skin  friction  is  independent  of  the  Mach  number 
and  depends  only  upon  the  surface  temperature . 

Heat- transfer  characteristics.-  The  heat-transfer  rate  per  unit 
area  at  the  surface  q  is  given  by  the  expression 


(-1) 


PrJvix 


M-  dfl 


and  from  equation  (l6) ,  since  3=1  and  9=0  for  tj  =  0, 


cPlTl  jp^Jh  7  -  1  2  *e  "  V  _1_  /plUlUl  fe 

“  ‘  Prn  V  x  2  V(co)  Pr-i  V  x 


where  7  is  the  ratio  of  the  heat  capacities  at  constant  pressure  and 
constant  volume . 
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As  can  be  seen  from  equation  (19),  the  heat-transfer  rate  is  pro¬ 
portional  to  the  difference  between  an  effective  enthalpy  function  oe 

and  the  surface  enthalpy  function  t3  s  (or  the  difference  between  the 
effective  enthalpy  ie  and  the  surface  enthalpy  is) .  In  contrast, 

the  heat-transfer  rate  is  approximately  proportional  to  the  difference 
between  the  recovery  temperature  (or  enthalpy)  and  the  surface  tempera¬ 
ture  (or  enthalpy)  when  the  Prandtl  number  is  taken  as  a  constant.  The 
effective  enthalpy  function  and  the  recovery  enthalpy  function  (the 
value  of  t3  at  the  surface  for  an  insulated  plate)  are  presented  in 
figure  3  as  a  function  of  the  stream  Mach  number .  Both  -0r  and  de 
decrease  with  increasing  Mach  number  and  de  decreases  with  increasing 
values  of  the  surface  temperature;  the  functions  i3r  and  i3e  become 
identical  for  the  condition  of  zero  heat  transfer. 


The  temperature  recovery  factor  r  may  be  determined  from  $r 

with  the  use  of  enthalpy  tables .  The  recovery  factor  is  presented  in 
figure  4  as  a  function  of  Mach  number  together  with  the  values  found 
from  other  analyses .  With  the  assumption  that  the  Prandtl  number  is 
independent  of  temperature,  the  recovery  factor  is  approximately  equal, 
to  the  square  root  of  the  Prandtl  number  for  laminar  flow  for  all  Mach 
numbers  (refs.  2  and  5,  for  example)  irrespective  of  the  viscosity- 
temperature  relation;  thus  the  recovery, factor  is  unity  for  a  Prandtl 
number  of  1.0  and  the  recovery  temperature  equals  the  stagnation  tempera¬ 
ture  of  the  stream.  When  the  variation  of  Prandtl  number  with  temperature 
is  taken  into  account,  however,  the  temperature  recovery  factor  decreases 
substantially  with  increasing  Mach  number .  The  recovery  factors  of 
reference  7 >  based  upon  the  use  of  the  Sutherland  equation  for  viscosity 
and  experimental  values  of  Prandtl  number,  are  in  substantial  agreement 
with  the  results  of  the  present  analysis ;  the  recovery  factor  is  approxi¬ 
mately  0.85  at  low  Mach  numbers  and  decreases  to  approximately  0.76  and 
0.65  at  Mach  numbers  of  5' and  10,  respectively. 

Equation  (19)  for  the  heat-transfer  rate  is  considerably  more  com¬ 
plicated  than  those  determined  from  analyses  based  upon  the  use  of  a 
constant  Prandtl  number.  The  function  de  is  dependent  upon  both  the 

surface  temperature  and  the  stream  Mach  number  for  a  given  stream  tempera¬ 
ture;  whereas  the  recovery  factor  (a  constant  for  all  Mach  numbers)  takes 
the  place  of  when  the  Prandtl  number  is  taken  as  a  constant . 

The  empirical  relation 


de  -  0S 


v(»)(dT.  -  o 


=  0.30532  -  0.001289 


2.25  -  M- 


1.35 


-  0.006062  Tg*  (20) 


i* 
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represents  the  variation  of  $e  with  within  1  percent  for  Mach 

numbers  up  to  5  for  a  stream  temperature  of  -67°  F .  Values  from  equa¬ 
tion  (20)  are  compared  with  the  calculated  values  of  the  quantities  in 

i9  _  ■d 

figure  5.  Eliminating  — - -  between  equations  (19)  and  (20)  leads 

V(“) 

to  the  following  expression  for  the  heat -transfer  rate: 

q  =  ?plTl  7  -.1  M2  (o.30532  -  0.001289 1 2. 25  -  m|1,35  - 

Pr±  V  x  2 

0.006062  Ts*  )(^r  "  ^s)  <21> 


The  values  of  $T  presented  in  figure  3  and  the  values  of  the  tempera¬ 
ture  recovery  factor  presented  in  figure  4  for  Mach  numbers  between  5 
and  10  were  determined  from  equation  (20)  with  values  of  i3e  and  V(°°) 

calculated  for  Tg*  =  2.  With  these  values  of  i3r  the  heat-transfer 
rate  calculated  from  equation  (21)  agrees  with  those  computed  from  equa¬ 
tion  (19)  for  Ts*  =  1,  2,  and  4  to  within  li  percent  for  Mach  numbers 

from  1  to  10 . 

Qomparison  of  heat-transfer  coefficients  or  Nusselt  numbers  found 
from  different  analyses  is  difficult  because  of  the  different  tempera¬ 
ture  functions  which  arise  in  the  expressions  for  the  heat-transfer  rate 
A  direct  comparison  of  the  heat-transfer  rates  more  clearly  shows  the 
effects  of  the  various  approximations  on  the  calculated  heat-transfer 
characteristics  •  The  heat-transfer  rate  determined  by  the  method  of 
this  analysis  is  presented  in  figures  6(a) ,  6(b),  and  6(c)  for  values 
of  Ts*  =  1,  2,  and  4,  respectively,  together  with  the  results  of  other 

analyses . 

Figure  6  shows  that  the  local  heat-transfer  rates  calculated  with 
the  use  of  the  Sutherland  equation  for  viscosity  and  a  constant  Prandtl 
number  of  0,75  are  in  very  good  agreement  with  the  results  of  this 
analysis  for  all  Mach  numbers  and  surface  temperatures  considered.  Like 
wise  the  results  of  reference  8,  which  employ  values  of  the  thermal 
properties  taken  in  part  from  experiment,  are  in  good  agreement  with 
the  results  of  this  analysis  for  Ts*  =  1,  the  only  value  of  surface 

temperature  where  a  comparison  could  be  made .  The  calculated  heat- 
transfer  rates  which  employ  more  approximate  viscosity-temperature 
relations  are  considerably  higher  than  the  values  determined  by  the 
method  of  this  report  at  very  high  supersonic  Mach  numbers;  below  a 
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Mach  number  of  4  or  5,  however,  the  agreement  is  generally  satisfactory 
for  a  stream  temperature  of  -67°  F . 


CONCLUDING  REMARKS 


Friction  and  heat-transfer  characteristics  of  the  laminar  boundary 
layer  on  a  flat  plate  have  been  calculated  for  Mach  numbers  from  1 
to  10  for  a  stream  temperature  of  -67°  F  with  the  use  of  experimental 
values  for  the  heat  capacity,  conductivity,  and  viscosity.  The  conse¬ 
quent  variability  of.  the  Prandtl  number  throughout  the  boundary  layer 
causes  the  recovery  factor  to  decrease  from  a  value  of  approximately 
0.85  at  low  Mach  numbers  to  approximately  O.76  and  0.65  at  Mach  numbers 
of  5  and  10,  respectively;  whereas  the  recovery  factor  is  approximately 
equal  to  the  square  root  of  the  Prandtl  number  for  all  Mach  numbers 
when  the  Prandtl  number  is  independent  of  temperature .  Moreover,  the 
heat-transfer  rate  is  proportional  to  the  difference  between  an  effective 
enthalpy,  which  is  a  function  of  both  the  surface  temperature  and  stream 
Mach  number,  and  the  surface  enthalpy  when  the  temperature  variation  of 
all  the  thermal  properties  is  considered;  whereas  the  heat-transfer  rate 
is  approximately  proportional  to  the  difference  between  the  recovery 
enthalpy  and  surface  enthalpy  for  calculations  based  upon  a  constant 
Prandtl  number .  Because  of  the  excellent  agreement  between  the  results 
of  the  present  analysis  and  calculations  based  upon  the  use  of  the 
Sutherland  equation  for  viscosity  and  a  Prandtl  number  of  0.J5,  it 
appears  that  the  exact  variation  of  the  thermal  properties  of  the  fluid 
need  not  be  considered  in  the  calculation  of  skin  friction  or  heat 
transfer.  More  approximate  viscosity-temperature  relations  and  the  use 
of  a  Prandtl  number  of  unity  introduce  significant  errors  in  the  calcu¬ 
lated  friction  and  heat-transfer  characteristics  at  Mach  numbers  above  1+ 
or  5  for  a  stream  temperature  of  -6 7°  F . 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  January  23,  1953 . 
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(a)  Insulated  surface . 

Figure  1--  Velocity  and  temperature  profiles. 
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M  Ts*  =  1. 


Figure  1.-  Continued. 
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p  Extrapolated,  from 
equation  (20)  for  Tg*=2 


Free- stream  Mach  number, 


Present  analysis 


-  Sutherland  eq.;  Pr=0.75  (ref, 5) 

-  p*=cT#;  P1-O.72  (ref. 2) 

- Extrapolated,  from  eq,(2C)  for  Tg*=2 

-  Sutherland  eq,,  Pr  variable  (ref. 7) 
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Free- stream  Mach  number, 


(a)  Tg*  =  1. 

Figure  6.-  Variation  of  local  heat-transfer  rate  with  stream  Mach  number. 
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(b)  T  *  =  2. 
x  7  s 

Figure  6.-  Continued. 
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